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COMPUTATION OF THE TWO-DIMENSIONAL FLOW IN A LAMINAR BOUNDARY 

LAYER 

By Hugh L. Dbtdbn 



SUMMARY 

A comparison is made of the boundary-layer flow com- 
puted by the approximate method developed by K. Pofd- 
hausen with the exact solutions which have been published 
for several special cases. A modification of Pohlhausen's 
method has been developed which extends the range of 
application at the expense of some decrease in the accuracy 
of the approximation. 

The work was carried out at the National Bureau of 
Standards, inpartwiththe cooperation and financial assist- 
ance of the National Advisory Committee for Aeronautics. 

INTRODUCTION 

The concept of the boundary layer introduced in 
modern aerodynamics' by Prandtl has been extra- 
ordinarily fruitful in the interpretation of experimental 
data. As yet, it is not possible to make the interpreta- 
tions quantitative, except in a few instances, since the 
equations describing the flow are nonlinear, and their 
mathematical solution is extraordinarily difficult, if not 
altogether impractical in many cases of interest. 
Pohlhausen (reference 1) developed an approximate 
method of solution of the equations for 2-dimensional 
laminar flow which has been criticised by von Mises 
(reference 2). Since Pohlhausen's method and related 
methods are within the mathematical skill of most 
experimenters, it seemed worth while to study the 
possibilities and limitations of such methods as judged 
by the instances for which exact solutions are known. 

POHLHAUSEN'S SOLUTION 

The equations for the steady laminar flow of an in- 
compressible fluid in the boundary layer along a 
2-dimensional surface whose radius of curvature is 
large as compared with the thickness of the layer are 
as follows: 1 
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where u is the tangential component of the velocity, 
v the normal component, x the distance measured along 
the surface, y the distance measured normal to the 
surface, v the kinematic viscosity of the fluid, p the 
density of the fluid, and p the pressure. At the surface, 
u and v are zero. As y increases, u approaches U, the 
speed in the potential flow outside the boundary layer, 
asymptotically. U is a function of a; in general. Prom 
equation (2) the pressure within the boundary layer is 
independent of y and equal to that in the potential 
flow. Since in the potential flow p + %pU* is constant, 



U' denoting^- 



(4) 



Prom equations (1) to (4), an important integral 
equation may be derived (reference 1), namely, the 
Karman integral relation 

-2U' Jj(U-u)dy- U± j\u-u)dy 

In Pohlhausen's approximate method of solution, a 
suitable assumption is made as to the shape of the 
velocity-distribution curve, leaving undetermined a 
parameter 5 which may be regarded as the "thickness" 
of the boundary layer. 5 is then determined as a 
function of x from the relation (5), following which the 
velocity distribution itself may be computed. The 
procedure is reviewed here, omitting the algebraic 
manipulations which are straightforward, though 
tedious. 

Pohlhausen assumes 

u^ay+btf+cys + dy* (6) 

To determine the 4 coefficients, 4 conditions are neces- 
sary. It is first required that the distribution within 
the boundary layer be continuous both as to magnitude 

' A fall discussion of the approximations made in deriving the boundary-layer 
equations is given by K. Hiemeni In Dlnglers Polyteohnlsohe Journal, vol. 326, 
911, p. 321. 
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and slope with, the potential flow at y — 5, where 5 is a 
function of x, at present undetermined. This gives the 
2 conditions: 

u= U at y=8 (7) 



bu n , 
^=0aty»5 



(8) 



The distribution is then made to satisfy the differen- 
tial equation (1) at the two boundaries. This 
requires * 

dy 5 pv Ox 

and 



at#=0 



(9) 



(10) 



From the 4 conditions 7, 8, 9, 10, the 4 coefficients 
a, b, c, din. (6), may be determined. It is found that 

~K)?' H-W; K—P 

d=(l — ^^r; that is, equation (6) is of the form 



(oc + o,X) f + (6 0 + &iX) ^ + (Co + c x X) ^ 



where X« 



Z7'5 2 



+ (d 0 +<Z 1 X)fJ 



(11) 



and do, ai, etc., are numbers having the 



values 6to=2, Oi = %, J 0 = 0, 6i=— K, — 2, d = K, 
do = l,c?i= — The parameter X gives the influence of 
the potential flow on the shape of the curve relating the 
nondimensional quantities u/U and y/d. If X is con- 
stant (or zero), the distribution curves are homologous 
for all values of x, i.e., the curve of u/U vs. y/6 is inde- 
pendent of x; S is, however, a function of a;. If 
do + Oi X is negative, u/U is negative for small values of 
y/6, indicating a reverse flow near the surface. The 



criterion for the beginning of reverse flow (separation) 
is that 

v Oi 



(12) 



The value of 5 is determined from (5), a procedure 
which amounts to satisfying the differential equation 
on the average and at the boundaries (by 9 and 10) but 
not at every pomt. "With the approximation (1 1) for the 
distribution within the layer, the \ipper limits of the 
integrals in (5) may be taken as 5 instead of <x> ; since 
at values of y/5>l, u is assumed equal to U and ?7-w 
vanishes. 

Trom (11) it may be shown that 



r 



(Z7— u)dy= US(S+ TX) 



(13) 



where S and T are numbers computed from a 0 , a it etc. 
Noting that both 5 and X are functions of x and that 

. U'S*. d\_2U'S dS , Z7"5 2 _2X dS , U"\ 
X= ~'dx" - ' 



dx 
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(14) 



where U" denotes 



^j o \u-u)d y =iS+3T»U$ 

+ T^-S\+U'6(S+T\) (15) 
Likewise JJ(C/-u) 2 dj/= U 2 8 (If+iX + MX 2 ) (16) 
whr re 2T,i,andMare numbers computed from Oo, Oi,etc, 
^ J\u-uydy= CP^ {K+ZL\ + bM\*) 



\U" 



udz -ao-QK-ZS+aJU'z- 
"2 dx = 



2i-3r+i 



+ 2UU'S(K+L\+M\ X )+U 2 S^(L+2MX) (17) 
Substitution in (5) gives, writing ^-z 

(£- T) ?7'V-2M(l +^P) U"z> 



K- S+Z (L- T) U'z+5MU' 

The values of S, T, E, L, and M are given by 
Of, bo Co do 



(18) 



T=- 



2 3 4 5 

o-i bi C\ di 

2" _ 3"~4 5 ~ 



■0.3 

1 
120 



■v i , «o 2 , V , Co 1 , do 2 „ 2&o Co 2do , cto&o 

20000 , gpgo , &qCq i 2Mo , Cpa\) 23 
+ 5" t "3' + "3" t "7" 1 "4 126 

r 20^ , 2&q6i , ?CqCi , "dpdi „ 26! d 2dx 

L ~— + -5~ + ~T + ~$ ai ~-3~~2~~5~ 

, gp&i I oA , 2aoC! 20,00 , gpfli , aidp 5oC t 
+ 2 2 + 5 5 3 3 3 



i Note that at jr— 0, u=r=0; at s-S, u- U, gj- ^.5^-0. 



,6lCo ,2Mi ,2M, Codj Cido 11 

+ !T + ~7 _ + 7 + 4 + 4 ~ 1512 



, , a? 61 2 Ci J d! 2 , a A , 2oiCi Oidj &iC, 

+ 7 + 4 9072 
The complete solution of Pohlhausen is given then by 

where X = u z = 

3 and hence 5 are to be determined from 
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dz 0.8 [-9072 + 1670.4X-(^ 



dz 
dx 



47.4+4.8 



UU" 
U" 



X 3 - 1 + 



UU" 

U' 



f7[-213.12 + 5.76X+X"J 



(20) 



Equation (20) is of the form 

dz P (x, z) 
3x = Q (x, z) 



(21) 



and cannot in general be directly integrated. In any 
particular case, a graphical solution can be made by 
the isocline method (reference 3) as illustrated later. 

MODIFIED POHLHAUSEN METHOD 

In the application of Pohlhausen's method to certain 
types of problems, difficulty arises because Q(x, z) 
vanishes and dz/dz becomes infinite. In at least one 
case in the literature, this behavior has been taken as 
an indication of early separation of the flow (reference 
4), but investigation shows that the singularity is 
introduced as a consequence of the nature of the 
assumed velocity distribution and represents a failure 
of the method. An attempt was therefore made to 
remove this limitation on the range of application of 
Pohlhausen's method. The attempt was not com- 
pletely successful; but, without serious additional 
complication, the range has been greatly extended, at 
the expense of some decrease in the accuracy of the 
approximation. 

The modification introduced was the addition of 
another term in the expression (11) for the velocity 
distribution, determining the additional constants to 
avoid the infinite value of dz/dx, if possible. Thus the 
velocity distribution was assumed to be 

g= (ao + a,X) | + (b 0 + WO ^ + (co + c,X) |J + (^ + d{K) ^ 

+ (eo+e,X)|J (22) 

Leaving Oo and di undetermined and applying the 
conditions 7, 8, 9, and 10, we find 

1 3 
6 0 = 0, h= co = 10-6ao, C! = 2-6oi, <4= -15 + 8<io, 



3 1 
di— — g + Sa!, 60 = 6 — 300, fii— ^ —3a!. 



Likewise: 



S-- 



l_Oo 
'2 10 



rp _1 <h 

120 10 

K= (2715 - 933oo + 104oo 2 )/6930 
L = (281 - 69oo - 3732oa + 832o o a 1 )/27720 
M =.(416o 1 2 - 69a! + 3)/27720 

The solution for X=0 (Z7' = 0) is known, namely, 
that given by Blasius (reference 5). The value of On 
was taken as 1.89 to give a good approximation to 
that solution. 

With tins value of a^: 
Q(x, z) =3328.4064-X(1837.44o a -241.23)-X 2 (2080ai I 
-3450a + 15) (23) 

To prevent dz/dx from becoming infinite, Q(x, z) 
must not vanish. Since in the physical problem X 
must be a real number, it would be desirable to have 
the roots of Q(x, z) = 0 imaginary. It proves to be im- 
possible to make the roots of #(a;,z) = 0 imaginary by 
any choice of <h. Calling the roots \ x and X 2 , where 

Xj >Xj and calling the value for separation xY = — 

<h was selected to make r 2 - as large as possible. 3 The 

maximum ratio is found for di = 0.11, X,= — 17.18, 
X 2 = -30.89, X! =48.52. The values for Pohlhausen's 
solution (21) are X,= -12, \ = -17.76, X 2 = 12. The 
range of application is thus extended by the modified 
method. The improvement in range is not accurately 
indicated by these figures since the values of 5 and 
hence of X are not strictly comparable; 4 but the im- 
provement is sufficient to deal with problems that 
cannot be handled by (21). 
The modified solution is then 

^= (1.89 + 0.11X)| - (|)fl+ (- 1.34 + 0.84X)fJ 

+ (0.12-0.62X)j£ + (0.33 + 0.17X)|£ (24) 



where X= U'z* 
mined from 



Z7'a 2 



z and hence 5 are to be deter- 



dz 
dx 



0.8[ - 59051 .9 + 13783 .3X - (53.93 - 14.6947 ^T)x 2 - (l + ^pr)* 3 ] 



U[- 1500.63 - 17.6337X + X 2 ] 



(25) 



» Tho modified method was developed in connection with the study of a flow in 
which separation was expected and In which Pohlhausen's solution failed by dj/dx 
becoming infinite. Hence this choice of m. The TrmTlmrrm ratio and the corres- 
ponding value of oi were found by trial, Le., by substituting various values of at In 
Qbc, t) <=0 and computing ^2. 



1 A more reliable Index of the improvement Is given by a consideration of the so- 
caHed "Verdrangungsdlcte", J*, defined by J*— J * M— jj \ d? and the corres- 
TJ'J*» 

ponding X*=— ^— • For Pohlhausen's solution, a*- 3(0300— O.00S33X), whence 

X*=X(0.300-0.00833A)'. For the modified solution S*-a(0.311-0.00287X), whence 
X*=X(0.811— O.O0267X) 1 . It is readily shown that for Pohlhausen's solution, X,* = 
-1.92, Xi'=-3.M, Xi*=0.<8; for the modified solution, X,*-— 2.19, Xi«=-4.78, 
Xi*-1.60. 
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APPLICATION TO PLOW IN WHICH U=kx- 
POHLHAUSEN'S APPROXIMATION 

A problem for which exact solutions have been given 

by Falkner and Skan (reference 6) is that in which 

U^hf where k and m are numerical constants. For 

this case TJ'=mTJIx, V =m(m-l) U/s?, and UU"/U'* 

= (m—l)/m, a constant. Equation (20) becomes, on 

dz 1 dX TJ" 
setting z=\/U', ^=-^7^-^75 X and collecting terms 

5 dX (X-XQ (X-X 2 ) (X-X3) , 2fi , 
3m+l X 3x~ X' + 5.76X- 213.12 ( OJ 

where X,, X 2 , and X 3 are the three roots of 

X 3 (0.6m + 0.2)+X 2 (36m + 1.92)+X(-213.12 

-1123. 20m) + 7257.6m = 0 (27) 

The variables X and x are separated, and equation (26) 
may be readily integrated to give * 

s • 

(X-X 1 )°(X-X 5 ) 6 (X-X 3 )' = ^ 5D,+I 
where A is the constant of integration, and a, b, and c 
are given by the relations 

X t * +5.76Xi -213. 12 
(X,-X s ) (Xi-X,) 

X 2 * + 5.76X 3 -213.12 
(Xs-Xa) (X s -X0 



a=- 



c= 



X 3 2 + 5.76X 3 - 213.12 
(Xs-Xj) (Xj-Xj) 



The constant of integration A must be determined 
by the boundary condition which fixes the value of X 
at some known value of 2. Thus if X=X C at x=l, 

A= (X c — Xi) a (X c — X 2 ) 6 (X e — X 3 ) c and the general solution 

fr-X.) a fr~*3) c „ ^ 

(Xc-XO" (Xc-X,) 6 (X^Xs)" 

The behavior of the general solution can be traced 
by somewhat tedious numerical calculations for definite 
numerical values of X e and m. If X e is selected equal 
to Xi, X 2 , or X3, it is obvious that the general solution 
degenerates to the particular solution X=Xi, X=X 2 , or 
X = X 3 , in which X is constant and independent of x. 
These are the so-called homologous solutions studied 

by Falkner and Skan. Smce X= = 

J V vx 



for this 



case, 5= -^/^^j The velocity distribution curves are 
homologous, the solution for any value of m being 
given by a curve of jj vs. V-^^^ 0T > s™ 106 ~ is constant 



for any m, of ^-vs. 2/-^/^." I* is convenient to write 
^■="jV, in which case equation (27) becomes 



* Assuming m^f—H, In which case there are only 2 roots of (27). 



iV 3 (0.6m 3 + 0.2m s )+iV 1 (1.92m+36m s )+iV(-213.12 
—1123.20m) + 7257.6 = 0 (28) 

Since (28) is a cubic in N, there are in general 3 
values of N and hence 3 solutions for each value of m. 
The results can best be visualized from a graph of N 
vs. m. This graph is shown approximately in figure 1 . 
The computed values are given in table I. 

TABLE I.— ROOTS OF EQUATION (28) 



m 






N, 








ZOO 


4.148 


8.088 


-38.64 


8.298 


16.170 


-77.28 


L60 


6.206 


1L200 


-6a 30 


7.808 


lasoo 


-76.46 


LOO 


7.052 


17.803 


-72.26 


7.062 


17.803 


-72.26 


aw 


1L210 


39.70 


-13a 60 


6.605 


19.860 


-06.25 


0 


34.054 


Infinite 


infinite 


0 


28.13 


-37.73 


-a os 


37.90 


— L46L46 


L 742. 70 


-a 768 


29.229 


-34.86 


-0.09 


75.28 


-386.57 


2ia84 


-6.776 


34.791 


-18.98 


-a 0929 


sa 68 


-378.1 


19L16 


-7.495 


36.13 


-17.76 


-a 10 


12a 00 


-m 0 


12a 00 


-12.000 


36.00 


-12.00 


-a 30 


Imag. 


-14,754.0 


Imag. 


Imag. 


442.62 


Imag. 


-0.33 


imag. 


Imag. 


Infinite 


imag. 


Imag. 


Inflnlto 


-a so 


Imag. 


Imag. 


362.3 


Imag. 


Imag. 


-181.16 


—LOO 


Imag. 


Imag. 


107.9 


imag. 


imag. 


-107.86 



•Xi-mJVi, X>=«niVi, Xj=mAV 

It may be noted that k does not appear in the solu- 
tion, k may be either positive or negative; equation 
(28) is the same in either case. If k is negative, U is 
negative, i.e., directed in the opposite direction to x 
and the negative values of N must be selected, since z 
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Fioube 1. — Oraph of equation (28), In part schematic 

is always a positive quantity, being proportional to 
the square of the thickness of the boundary layer. 

The solution given by Pohlhausen for the case m = 
— l,k negative, is wrong. As shown by figure 1, 
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there is no solution given for this case by Pohlhausen's 
method. The solution given in reference 1 was appar- 



ently obtained by setting 



— x s g^ + 2Xa;, which 



is obviously wrong. 

For the flow i7=fcc", no difficulty is encountered be- 
cause of the vanishing of Q(x,z) in equation (21). X is 
constant for a given value of m and when m is such 
that the corresponding \( = mN) is equal to a root of 

dz 

Q(x,z)=0, X is also a root of P(z,z)=0. Hence ^ 

is indeterminate, but not infinite. However, when X 
is greater than 12 (one of the roots of Q(a;,z)'= 0) the 
speed within the boundary layer rises to a maximum 
exceeding U and then falls to U. Although such solu- 



to the assumptions on which the approximate equations 
were deduced. 

This leaves for consideration the branch labeled A, 
divided in three sections A u A 2 , and A 3 . Ai and A 2 
represent flows without separation, X being between 
+ 12 and — 12. The branch A3 represents a flow with 
separation. For negative values of m and positive k, 
representing a flow with velocity decreasing as a; in- 
creases, if m does not exceed in absolute value 0.10, 
there are 2 solutions, 1 without separation and 1 with 
separation. For larger negative values of m, there is 
either no solution at all or only a solution with separa- 
tion. This result may be compared with that of Falk- 
ner and Skan where a flow without separation was 
possible if m was negative but not for absolute values 




— Folhner and Shan 
—Pohlhousen 

*7 Modified PohlhauserT 

+ Blosius 



3 4 5 6 

yT/U/xv 

Figure 2.— Velocity distribution In boundary layer for the case UMfcr- 



tions seem to be possible even in the exact treatment 
given by Falkner and Skan, they do not, so far as 
known, occur in any actual flow. The solutions repre- 
sented by the branches marked B in figure 1 (N 3 and 
Xj in table I) are therefore probably not of interest, 
likewise when X is negative and much greater in abso- 
lute value than the other root, — 17.76, large negative 
values of u occur within the boundary layer. The 
solutions represented by the branches marked O in 
figure 1 (N 3 and X 3 in table I) are probably not found in 
any actual flow. It must be remembered that although 
such solutions of the boundary-layer equations- may 
exist, the boundary-layer equations are themselves 
approximations. The solutions represented by 
branches B and O are of the type which do not conform 



greater than 0.09. The branch A lt represents a flow 
with a velocity increasing as x increases. 

The speed distributions for positive values of k and 
values of m equal to 2.0, 1.0, 0.5, 0, and —0.09 were 
given by Falkner and Skan. They are reproduced in 
figure 2, together with the results computed from table 
I (root Ni) and those obtained by the method given 
later. 

The case m = 0 is that treated by Blasius, whose 
results are also shown. 

The agreement is very close except for m=» —0.09. 
For this case the method of Falkner and Skan is prob- 
ably open to criticism. The series used to represent 

the solution is not convergent for values of y-^^mach 
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greater than 6, at which, value u/U equals 0.95. At 
y-^^=5, u/U is 0.83 and it appears impossible to tell 

whether u/U approaches 1 as y increases, since the 
series is not convergent. Since the approach of u/U 
to 1 is the criterion for determining the constant which 
deterrnines the coefficients in the series, it cannot be 
demonstrated that the solution given is correct. How- 
ever, this difference may be taken as a warning that 
Pohlhausen's method may not be satisfactory for 
negative values of U'. a 
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Fiqube 3. — Graph of equation (29), in part schematic 
MODIFIED METHOD 

The solution by equation (25) proceeds along the 
same lines as by Pohlhausen's approximation. Equa- 
tion (28) is replaced by 

iV 3 (0.6m 3 +0.2m s )+iV i (49.02194m 2 -5.87794m) 
+iV(-9526.01m-1500.63) +47241.52 = 0 (29) 

The roots are shown graphically in figure 3 ; the data 
from which the graph is plotted are given in table II. 



• Some recent experimental work completed at the National Bureau of Standards 
Indicates that Pohlhausen's approximation Is indeed very poor for negative values 
of V, separation occurring at values of X of the order of — S. 



Tablb H.— ROOTS OF EQUATION (29) 



n 


Ni 


Nt 




Xi 


Xi 


Xi 


2.0 


2.352 


44.797 


-8a066 


4.704 


8a 69 


-16a 13 


L5 


3.056 


ea 009 


-104.067 


4.684 


9a oi 


-166. 10 


L0 


4.366 


9a 76 


-149.06 


4.366 


9a 70 


-149.00 


as 


7.638 


186.12 


-267.29 


3.819 


92.66 


-133.66 


0 


3L481 


Infinite 


Infinite 


0 


102.66 


-73. 16 


-a ofl 


84.956 


377.61 


-1,246.6 


-7.646 


11410 


-33.88 


-a 10 


117.61 


251.93 


-1, 139.4 


-11.761 


113.94 


-28.10 


-a 1036 


169.8 


169.8 


-1,106.9 


-17.693 


114.69 


-17.69 


-a is 


Imag. 


Imag. 


-839.3 


Imag. 


Imag. 


164.28 


-a 20 


Imag. 


Imag. 


-862.9 


Imag. 


Imag. 


17a 68 


-a 30 


-43.66 


-177.42 


-36,086.0 


13.066 


63.23 


10,620.0 
Infinite 


-0.33 


-33.03 


-193.09 


Infinite 


1L010 


64.30 


-a so 


-16.661 


-166.086 


778.6 


7.826 


77.64 


-388.3 


-LOO 


-6.168 


-84.26 


227.7 


a 168 


84.25 


-227.7 



The branches A u At, A3, B, and C correspond to 
similar branches in Pohlhausen's solution. The new 
feature is the occurrence of branches A4 and Bu B\ 
is of the B type discussed previously. Ai gives a solu- 
tion without separation for negative values of m greater 
in absolute value than 0.276, if k is negative; i.e., if the 
flow is one in which the speed increases as x increases. 
For m=— 1, the boundary-layer equation may be 
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Fioubb 4.— Veloolty distribution In boundary layer for the case U=-*/r 

exactly integrated (reference 1) with the resultant dis- 
tribution: 

9f-2-8tanh«(l.l46+i^3) 

It should be noted that since k is negative, U is nega- 
tive, and hence — U is positive. Whereas equation 
(28) gave no solution, equation (29) gives 



US' 

N-= — = -6.158, X=6.158, 3 = 
vx 



= 2.481 



Vvx 



The corresponding velocity distribution from (24) is 
compared with the exact distribution in figure 4. 
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The results for m equal to 2.0, 1.0, 0.5, 0, and -0.09 
are shown in figure 2 for comparison with the results by 
Pohlhausen's solution (equation 28) and the results of 
Falkner and Skan. It is seen that Pohlhausen's solu- 
tion is in general a better approximation than the 
modified solution. 

In figure 1 and figure 3 the values of X change con- 
tinuously along the several branches of the curves. 
The way in which the branches join at infinity is indi- 
cated by the variation of X. A study of these two figures 
leads one to believe that equations such as (20) and 
(25) of the third degree in X cannot represent the solu- 
tion over the entire range, and since the degree in X 
can be traced to the fact that (19) and (24) are linear 
in X, it is probable that the distribution curves either 
axe not linear in X or that other quantities such as 
UU"\U' % also affect the shape of the distribution curve. 
Nevertheless, the approximations are valuable where 
they do give a solution. 

APPLICATION TO TRANSVERSE FLOW ABOUT A 
CYLINDER 

The flow in the boundary layer of a cylinder has been 
computed by J. J. Green (reference 7) by a step-by-step 
method in which it was assumed that the circumferen- 
tial velocity is expressible as a power series whose 
coefficients are functions of the distance along the sur- 



speed U 0 , x the ratio of the distance measured along 
the surface from the front stagnation point as origin 
.to D, y the ratio of the distance measured normal to 
the surface to D and 5 the ratio, of the thickness of 
the boundary layer to D. It may be shown that the 
only change required in (25) is in the definition of z 
which becomes z=*RS 2 where R is the Reynolds 
Number UJ)/v; in other words v is replaced by 1/R. 

The values of U, U' , and U" were taken from data 
given in Green's paper, table ILT, columns 4, 5, and 
8. The relations between U, U' , and TJ" and Green's 
fi,ji, and gR are as follows: 
E7=2 

U R'q 
UU" .RW 



U' 1 



— 1 



The solution of equation (25) was carried out by 
the isocline method. The first step in this method is 

dz 

the computation of values of ^ for several values of 

UU" 

z and x. The values of x, U, U', and -jj, s are given 

.dz, 



in table HI, together with the values of ^ for several 
values of z obtained by substitution in equation (25). 
Table EH. — DATA FOR SOLUTION OF EQUATION (25) FOR CYLINDER 



8 


X 


U 


W 


VTJ" 




0 


a 76 


1.50 


2.28 


3.00 


3.76 


o 

0 


0 


0 


a 892 


-a 2636 


+» 


+« 


— CD 


— 00 




— 00 


6 


aoiso 


a 149 


3.636 


+.0446 


212.0 


7a 6 


-49.4 


-169 0 


-288°0 




10 


.0878 


.307 


1816 


—.0270 


102.7 


34.0 


-29.86 


-9tt 4 


-isao 


"-mo 


16 


.1303 


.486 


3.643 


-.0681 


66.0 
499 


23.4 


-16.30 


-62.0 






20 


.1744 


.633 


3.636 


-.2746 


18.81 


-9.90 


-37.3 


-64.0 


-99 6 


28 


.2180 


.789 


3.243 


-.4236 


39.9 


17.07 


-4.21 


-24.6 


-44.2 




30 


.2616 


.928 


a 018 


-.4816 


33.9 


16.83 


-1.13 


-17.29 


-33.0 


^4sT~ 


36 


.3062 


1.066 


2.846 


-.4486 


29.8 


14.79 


aw 


-12.81 


-26.88 




40 


.3488 


1.176 


2699 


-.6387 


26.79 


13.97 


1.89 


-9 63 


-2tt 79 


-3L8 


45 


.3924 


1.283 


2.602 


-L4204 


24.64 


13.64 


3.36 


-a 47 


-16.96 




60 


.4360 


L384 


2,113 


-3.9477 


22.76 


14.21 


a 14 


-L68 


-9.03 


-1134 


66 


.4796 


1.463 


L643 


-a 1623 


21.61 


16.89 


9.96 


4.68 


-.62 


-6.67 


60 


.8232 


1.623 


.967 


-19 468 


2a 69 


17.16 


1177 


19 68 


7.44 


4.44 


66 


.6668 


1.649 


.462 


-9L04 


20.33 


18.69 


17.19 


16.78 


14.48 


13.26 


67 


.6838 


1.664 


.161 


-923.0 


2a 27 


19.73 


1984 


19 08 


1193 


18.93 


69 


.6016 


1.654 


-.096 


-3,748.0 


2a 27 


25.91 


26.73 


27.81 


29.16 


3a 77 



face. The approximate solution by equations (24) and 
(25) was computed for comparison with Green's more 
nearly exact solution. 

While the quantities ?7~ij~> ^> jj/i hi equa- 

tion (25) are nondimensional, z and x are not. As 
noted by Green, it is convenient to measure all dis- 
tances in terms of some reference distance, in this 
case the diameter D of the cylinder, and all speeds in 
terms of some reference speed, in this case the speed 
U Q at a great distance from the cylinder. For sim- 
plicity, no new symbols will be introduced, but U is 
taken to mean the ratio of the speed in the potential 
flow outside the boundary layer -to the reference 



In addition to x, which is the distance along the sur- 
face from the front stagnation point, the azimuthal 
angle 8 is given. 

An isocline diagram is prepared from the data in 
table m, that is, a chart with z as ordinate, and x as 
abscissa with curves showing the loci of constant values 

dz 

of • In practice it is convenient to change the scale 

dz 

of z relative to a; to give values of less than 10. 

In the present case the use of 'a/10 and^ is found 

desirable. A portion of the isocline diagram is shown 
in figure 5. The numbers on the curves are the values 
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The curves are located by interpolation be- 
tween the values given in table HE. Thus at x = 0.4796- 



the isocline 



dz/10 
dx = 



/dz 



') 



10 ) lies at a value of z of 



1.49; the isodine- ^^ « 



• 1.5 lies at 2=0.829, etc. 



The solution curves of equation (25-) must cross the 
isoclines with the slope indicated on the isocline, that 
is, the zero isocline must be crossed horizontally, the 
isocline labeled 1 at a slope of 45°, etc., as indicated by 
the short hues in figure 5 crossing the isoclines. The 
particular solution curve in which we are interested is 
the one which satisfies the boundary conditions at 
£■=0, the front stagnation point. Because U=0 at 
x= 0, we find here a singularity; there is no true bound- 
ary layer right at the stagnation point. The isocline 

dz 

chart shows a singular point at which^is indetermi- 
nate. No matter what value of z is assumed at x-=0, 



0.4 




0.3 0.4 0.5 
x 

Fiqube 5. — Isocline diagram for boundary layer of a cylinder 

the attempt to construct the solution curve leads imme- 
diately to the singular point. 

It may be shown that the zero isocline leaves the 
singular point at zero slope and that the desired solu- 
tion is constructed by starting a solution curve in this 
manner. The curve is shown in figure 5. In the actual 
computation, a greater number of isoclines were drawn 
to a larger scale. From this curve, the values of 
2= .55 \ were, obtained as. a function of x (table IV). 

Table IV.— VALUES OF Ri 5 FOR CYLINDER 



X 


Ri 1 


X 




0 


1.126 


0.3488 


1.600 


0.0436 


1.125 


.3924 


L690 


.0872 


L140 


.4360 


1.731 


.1303 


1.186 


.4796 


1.966 


.17+4 


1.233 


.6232 


2.340 


.2180 


1.300 


.6668 


2.880 


.2816 


1.377 


. .6838 


3.166 


.3052 


1.433 


.6016 


3.600 



The velocity distribution was then determined from 
(24). The curves for x= 0.1308 and 1=0.5668 only 
are shown in figure 6. The results of experiments by 



Green and the results of Green's computation are also 
shown. At other values of x, the differences are of a 
similar nature. In Green's computation, the speed in 
the potential flow was taken from the experimentally 
measured values, although the pressure distribution 
was also used in the remainder of the computation. Wo 
have used only the results of the pressure distribution. 

REMARKS ON THE ACCURACY OF THE APPROXIMATE 
METHODS 

The preceding comparisons show that in theso 
particular cases approximate methods give a fairly 
good representation of the actual distribution, the 
differences not exceeding 0.05 U for the modified 
solution, or 0.02 £7 for Pohlhausen's solution, where it is 
applicable. Unfortunately, all of the satisfactory 
exact solutions are cases in which X is positive and 
less than 10 and UU"/U' 1 is small. No satisfactory 
comparisons are known in which X is negative. 



1.0 



o.a 



u 



as 



0.4 



OB 





























> 


■ - c 


use 


8 * 
















— o — 










y 
























t 

t — 




'/ x 


= O. 


566t 


i 














II 
f * 


r 
























r 


























~fr — 
It 

ft i 
! £ 


























' ^ 








-een 
<-om 
■een 
odffl 






















(f 

12, 


a uieor&Ticai 
small scale plotj 












— o 
+ 




— M 


o < 

3d F 


to/?//- 


imat 

oust 


3D S 


oluh 


on 































0.4 



0.8 



2.0 



24 



1.2 1.6 

Fioube 6.— Velocity distribution In boundary layer of cylinder at two distances 
from the front stagnation point 

It is possible to approach the question in a different 
way. The approximate solution may be used to com- 
pute all quantities entering in the differential equation 
(1) and a check made as to the accuracy with which 
the equation is satisfied. There will be found a 
residual error which is most conveniently expressed as 

a ratio to the last term UU' by equation 4). 

This error will b.e a function of y, nearly zero on the 
average, exactly zero at y = 0 and y = 5 and at some 
intermediate point. Figure 7 shows the maximum 
positive and negative residuals for values of X from 
TJTJ" 

0, +44, and -45 for the 



-20 to +20 for 



UU' 



modified solution (equations 24 and 25). 

The maximum errors for X between 0 and 8, jjrr = 0, 
the region in which comparisons with exact solu- 
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tions are available, are of the order of 0.3 to 0.4 
times TJU'. (Note that X = 0 corresponds to Z7'=0, 
hence error jUU' approaches infinity at \ = 0). The 
errors for negative X are slightly greater than for 
positive X. As UU"/U' 2 increases, the errors become 
larger, especially for positive values of X. 
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Fiqube 7.— Errors In the modified solution (equations (24) and (26)) at various 
values of X and UXT'IIT' 

CONCLUSION 

Approximate methods of the type suggested by 
Pohlhausen for calculating the flow in a laminar 
boundary layer are useful in giving one a fair picture 
of the flow when the parameters X and UU"/U' 8 are 
not too large. The solution given by Pohlhausen fails 
when U' is positive and large, such that X = 12. An 
extension of the range of application of the solution has 
been accomplished by a modification of Pohlhausen's 



method with a decrease in the accuracy of the approxi- 
mation. Comparisons have been made of the approxi- 
mate solutions -with exact solutions for the cases in 
■which exact solutions have been published. The 
approximate solutions of the type studied appear to 
be very poor when X is negative. 
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